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Abstract 

We propose a new test for distinguishing observationally cosmolog- 
ical models based on seed-like primordial perturbations (like cosmic 
strings or textures), from models based on Gaussian fluctuations. We 
investigate analytically the Fourier space statistical properties of tem- 
perature or density fluctuation patterns generated by seed-like objects 
and compare these properties with those of Gaussian fluctuations gen- 
erated during inflation. We show that the proposed statistical test can 
easily identify temperature fluctuations produced by a superposition 
of a small number of seeds per horizon scale for any observational an- 
gular resolution and any seed geometry. However, due to the Central 
Limit Theorem, the distinction becomes more difficult as the number 
of seeds in the fluctuation pattern increases. 
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1 Introduction 



One of the directly measurable features of the primordial fluctuations that 
gave rise to galaxies and large scale structure formation is the probability 
distribution and the corresponding moments of the primordial perturbation 



The phases 4>k of the Fourier modes S(k) are ususally assumed to be mi- 
correlated and randomly distributed according to a uniform probability dis- 
tribution. This assumption, based on the prediction of inflationary models, 
leads by the Central Limit Theorem to a Gaussian probability distribution 
for the field S(x). 

There are two main advantages of such Gaussian models: First, all the 
statistical information about the field S(x) is encoded in a single function: the 
two point correlation function (or equivalently the power spectrum). Second, 
when combined with Cold Dark Matter (CDM), Gaussian models are in 
reasonable agreement with small and intermediate scale observations (White 
et. al. 1987). However, observations on large scales (larger than 10/i~ 1 Mpc) 
have consistently indicated that Gaussian CDM models lack power on large 
scales. 

One approach to the resolution of this problem is to retain the Gaussian 
nature of the primordial perturbations while modifying other aspects of the 
model in an effort to transfer power to large scales. This has led to the 
construction of the 'hydrid models' which attempt through the introduction 
of additional parameters (like a component of Hot Dark Matter (HDM)) to 
reconcile Gaussian models with large scale structure observations. 

The other approach, is the consideration of non-Gaussian primordial per- 
turbations. A class of non-Gaussian perturbations which is well motivated 
physically is seed-like perturbations. These primordial perturbations may be 
naturally provided by topological defects (e.g. cosmic strings (Kibble 1976; 
Vilenkin 1981; Brandenberger 1992) or textures (Turok 1989)) produced dur- 
ing phase transitions in the early universe. Other interesting mechanisms 
(e.g. primordial black holes (Carr & Rees 1984)) can also produce seed- 
like perturbations. Models based on cosmic strings for example, have been 
shown to have several interesting features that make them worth of further 
investigation. After appropriately normalizing the single free parameter of 
the model, cosmic strings can naturally provide concentrations of galaxies on 
sheets (string wakes) with typical dimensions 40 x 40 x Ah^Mpc 3 (Vachas- 
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pati 1986; Stebbins et. al. 1987; Perivolaropoulos, Brandenberger & Stebbins 
1990; Vollick 1992; Hara & Miyoshi 1993), they are consistent with the re- 
cent detection of anisotropy by COBE (Bouchet, Bennett & Stebbins 1988; 
Bennett, Stebbins & Bouchet 1992; Perivolaropoulos 1993a), and they are 
in reasonable agreement with observations of peculiar velocities regarding 
measurements of the Cosmic Mach Number (Perivolaropoulos & Vachaspati 
1993). However, like the CDM model, the cosmic string model is not free 
from problems. As pointed out by Albrecht & Stebbins (1992a) the power 
spectrum of density fluctuations produced by cosmic strings in a universe con- 
sisting mostly of CDM appears to have too much power on small scales. This 
problem was shown to be resolved however, if CDM is substituted by HDM 
(Albrecht & Stebbins 1992b). In addition, Perivolaropoulos & Vachaspati 
(1993) have recently pointed out that cosmic strings can not explain the ob- 
served magnitude of peculiar velocity flows on scales larger than 50/i _1 Mpc 
if normalized from peculiar velocity observations on scales 5 — 20h~ 1 Mpc. 
This problem however, also appears in the CDM model and may be resolved 
by assuming velocity bias. 

Perturbations in seed-based models may be represented as a superpo- 
sition of localized fluctuations with geometry that depends on the model 
under consideration. The most sensitive way to distinguish observationally 
models based on Gaussian perturbations from models based on seeds is by 
investigating the statistical properties of the perturbations. In fact, inter- 
esting statistical tests have been proposed that attempt to provide ways to 
efficiently make this distinction (Coles 1988; Lucchin, Matarrese & Vittorio 
1988; Scherrer, Melott & Shandarin 1991; Gaztanaga & Yokoyama 1992; Luo 
& Schramm 1992; Perivolaropoulos 1993b). 

However, there are two main problems that tend to decrease the sensi- 
tivity of these tests. The first comes from the Central Limit Theorem which 
predicts that as the number of superimposed seeds increases, the resulting 
perturbations look more like Gaussian. The second comes from the finite res- 
olution of observational experiments. Observations effectively average over 
patches in the sky and associate with each patch a measurement that may 
correspond to a temperature, a density or a velocity field. By the Central 
Limit Theorem such averaging tends to reduce the non-Gaussian signature 
of seed based models. 

The statistical test we discuss in this paper is an attempt to evade the 
second problem. By studying the statistical properties of Fourier modes we 
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can isolate the effects of low resolution, to high wavenumber k Fourier modes 
and focus on the low k modes that remain unaffected by the smoothing on 
small scales. 

In what follows we consider perturbation patterns produced by a random 
superposition of N identical seed perturbations and derive the probability 
distribution and moment generating function of the Fourier modes that cor- 
respond to the pattern. The pattern of perturbations investigated here is 
known in the literature as 'shot noise' (Campbell 1909; Rice 1944) and ap- 
pears in several and diverse problems. The statistical properties of shot noise 
have been studied previously (Rice 1944) mainly in coordinate space and in 
the large N limit, showing strong Gaussian behavior. In the present analysis 
we focus instead on the statistical properties in Fourier space. The results 
presented here are fairly general in that they are valid for any value of N 
and any shape of the superimposed seeds. For simplicity we focus on the one 
dimensional case but we show that the analysis can be easily generalized to 
higher dimensions. 

2 The Large N Limit 

Consider the random function 

f(x) = Y,f 1 (x-x n ) (1) 

n=l 

where fi(x) is a seed function superimposed randomly at positions x n such 
that —l<x n <l. Both f(x) and fi(x) are defined within the interval [— /, I] 
and periodic boundary conditions are used during the superposition. The 
Fourier expansion of fi(x) is: 

= E 9i(k)e l ^ x (2) 

fc=— oo 

with 

9i(k) = ± £ dxfMe-*** (3) 
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The fact that fi(x) is real implies that g\{k) = gi(—k). Using (1) and (2) 
the random function f(x) can be expanded as: 



/(*) = £ 9i(k)e^Q(k) (4) 

fc=— oo 

with 



N 



Q{k) = J2 e 1 1 Xn = <?iO; X!..x N ) + iq 2 (k; x x ..x N ) (5) 

n=l 

Thus all the random properties of f(x) have been transferred to Q(k) which 
is independent of the shape of the seed function fi(x) and can be viewed as 
the final position of a N step random walk in the two dimensional qi — q 2 
plane. Notice that the reality condition Q*(k) = Q(—k) which is trivially 
satisfied implies that the end points of random walks with negative k are 
simply obtained by reflection with respect to the q\ axis of the corresponding 
positive k end points. We are interested in the joint probability distribution 
P(qi,q2) and the corresponding moment generating function. 

The general case of arbitrary N is treated in the next section. Here we 
study the special case N — > oo for which there are results available in the 
literature (Rice 1944). For iV — > oo it is easy to show that the random 
variables qi,q 2 become independent i.e. 

P(qi,q2)=P(qi)P(q2) (6) 

In the same limit, the Central Limit Theorem implies that both qi and q 2 
(being sums of identically distributed random variables) are distributed ac- 
cording to the Gaussian 

1 fe-M) 2 

where \i =< qi > and a 2 =< qf > {i — 1,2). Since the probability distribu- 
tion of x n is uniform in the interval [— /, /] we can find \i and a 2 as: 

1 r+ l N kir 

/i =< qi >= (—) N J i dxi...dx N (J2 cos — x n ) = N5 k0 (8) 

and 

° 2 =< q\ >= y (9) 
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Similar results are also easily shown to hold for q 2 . Thus for k ^ we find 

P(gi,g 2 )->^e-£ (10) 

where q 2 = qf + q 2 - Clearly, the probability distribution is independent of 
the phase 0& = tan~ 1 ( 2 j) of the Fourier modes. Therefore, for N — > oo the 
Fourier phases are distributed uniformly while the Fourier mode magnitude 
q(k) has a Gaussian distribution. It may be easily seen by visualizing the 
random walk Q{k) that the probability distribution of (pk will in fact be 
uniform for any N. However, the rest of the results of this section are not 
valid for finite N since the independence of the variables 51,52 (expressed 
through (6)) breaks down in that case. This will be shown rigorously in the 
following section. 



3 Arbitrary N 

The Fourier transform of P(5i, 52) may be written for any N as: 

/+N 
d qi dq 2 P( qi ,q 2 )e^e 1 ^ (11) 
-N 

which implies that 

1 +00 

P(5i,5 2 ) = (^) 2 E P(Pi,P2)e-^e-^ (12) 

P1,P2=— OO 

where pi,p 2 are integer variables, Fourier conjugate to 51,52- By inspection 
of (11) it becomes clear that P(pi,p 2 ) is also the moment generating function 
for the distribution P(5i, 52). In fact, it is easy to see that 

< ft 12 >~ (^) dpTdp n l Pl = P2 =0 (13) 

But the same moments are also generated by the function 

R( PuP2 ) = (1-) N / +i rfa; 1 ...rfx J ve l ^ 9l(A;;:Cl '- :rjv) e i£ ^ 92(fc;a;i - a;iv) (14) 
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since h is the probability that x n will be in the range [x n ,x n + dx n ]. It 
is easy to check that (13) also holds with P(pi,p 2 ) substituted by R(p\,p 2 ). 
Since the moment generating function that corresponds to a given probability 
distribution is uniquely defined (Feller 1971) we must have 

P(pi,P 2 )=R( P l,P 2 ) (15) 

By expanding qi(k; Xi...x^) and q 2 (k; Xi...Xn) according to (5) we obtain 
using (14) and (15) 

P( Pl , P2 ) = (-L ^ e *(*icos£+t 2 sinO)iV ( 16 ) 

2-KK J-kw 

where U = p^ (i = 1,2) and £ = Let now 

h = tcos5 (17) 
t 2 = tsin5 (18) 

Using the periodicity of cos£ and the fact that k is integer (16) becomes 

P(Pi,P2) = (^J^d^ tco ^) N (19) 

or 

P(pi,P2) = Mt)f (20) 

where 

t = (jj)y/A+ti = (jj)p (21) 

The generating function (20) is one of the central results of this paper. It is 
valid for any N and clearly depends only on the magnitude p of the vector 
ipiiPi)- This implies that its Fourier transform P(qi,q2) will similarly be 
a function of the magnitude q only and there will be no dependence on 
the phase of the vector (qi,q 2 ). Thus, the Fourier phases of seed induced 
perturbations obey a uniform distribution for any number N of superimposed 
seeds. Obviously, this statement applies to each mode k individually and does 
not imply that there will be no correlations among the phases of different 
modes. Such correlations will clearly exist for seed perturbations but are not 
the subject of the present study. 
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Our results can explain the numerical simulations of Suginohara & Suto 
(1991) where it was found that even in strongly non-Gaussian evolved density 
fields the phases 0& are uniformly distributed. The authors of that paper had 
concluded that the investigation of the distribution function of the phases fa 
does not provide a sensitive test of the non-Gaussian behavior in the strongly 
non-linear regime but no clear explanation was given of this fact. Since the 
density field in the non-linear regime can be viewed as a superposition of 
dense lumps (seeds), the above analysis is applicable and predicts exactly 
the uniform distribution of phases seen in the simulations of Suginohara & 
Suto (1991). 

The probability distribution P(qi, 92) = P(q) is obtained by Fourier trans- 
forming the generating function (20) as follows: 

1 +°° 

n<7) = (^) 2 E (Jo(t(pi,P2))) N ^ (22) 

P1,P2=— OO 

where t = ^(pi,^) an d q = (91,92) ■ For N > 1 we may approximate the 
sum (22) by an integral and reduce it to 

! r°° 

P(q) = 7T dpp(J (p)) N J (pq) (23) 

27T JO 

It is straightforward to verify that for N > 1 

/oo roo roo roo 

d Ql / dq 2 P( qi ,q 2 ) = / dpp(J (p)) N / dqqUpq) = 1 (24) 
-00 J —00 jo Jo 

which is to be expected since P(9i, 92) is a probability distribution. 

Directly measurable quantities from a given fluctuation pattern are the 
moments of the fluctuation probability distribution. The moments generated 
by the function (20) can easily be obtained and compared with the moments 
of the Gaussian probability distribution. Since P{j)\,P2) depends only on the 
magnitude p, it is easy to show using (13) that for any positive integers m 
and N we have 

< 9i >-< 9 2 >- \t=o l^oj 

By expanding (J (t/i)) N in powers of t we obtain 

w/i)f = 1 + {t ^ y + {t4m)\\ - ^) + - (26) 
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From (25) and (26) it may be shown that the moments of the appropriately 
normalized variables tv = — % — (i — 1, 2) for k 4^ are 



< rf m+1 > = (27) 
<rf> = 1 (28) 

<r\> = 3(1-^) (29) 

The kurtosis (defined as (< r| > —3)) is negative for all finite N and ap- 
proaches the Gaussian value for large N. Also, the skewness < rf > is for 
all N. By expanding the generating function further, the higher moments may 
also be obtained. The negative sign of the kurtosis is to be contrasted with 
the corresponding sign of the kurtosis of seed perturbations in coordinate 
space where several cases of interest have been shown to have positive kur- 
tosis (Scherrer & Bertschinger 1991; Luo & Schramm 1992; Perivolaropoulos 
1993b). 

By numerically evaluating the integral (23) we plot the probability dis- 
tribution P(q) and compare it with the Gaussian. This is shown in Figure 1 
(dotted line) for N = 10. The corresponding Gaussian distribution with the 
same standard deviation (obtained from (10) with N = 10) is also shown for 
comparison (continous line). 

It is of interest to obtain the generating function for the moments of the 
normalized variables rj. This is easily shown to be 

P(h, h) = (M^—)) N = (i + ^7 + ...)" — ^ (30) 

iy/N/2 2N 

where the limit, indicated by the arrow, leading to the standard Gaussian 
generating function e 2 , is obtained for iV ^> 1. Thus, the generating function 
approaches, as expected, the Gaussian for large N. 

Let us demonstrate the utility of these results in a somewhat realistic 
case. Consider an one dimensional pattern of fluctuations in Fourier space. 
In a realistic case, these fluctuations will be a superposition of a Gaussian 
noise random variable q n and a signal q s . Let the signal to noise ratio be 
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assumed known. The measured variable at each pixel is q = q n + q s . We 
want to test the hypothesis that q s is produced by a superposition on N 
seeds. The moment generating functions for the variables q n and q s are: 



*2 



M qn (t)=e a ^ (32) 

and 



M q .(t) = (J (t/i)) N (33) 



Therefore 7 = < ^| > = Since the variables q n and q s are independent, the 
moment generating function for the measured variable q is: 

M qs (t) = (Mt/i)) N e^ (34) 

It is now straightforward to expand M q (t) and thus obtain the kurtosis k for 
the random variable q: 

K = T7^ = 3{1 -2(i + 1 yN ) (35) 

The kurtosis k is measurable, and any constraint on it can be translated using 
(35) to a constraint on N, the number of superimposed seeds on the pattern 
under consideration. Given that different seed-based models predict widely 
different number of seeds per Hubble scale (according to simulations, there 
are 0.04 textures unwinding per Hubble volume per Hubble time while the 
corresponding number for long strings is about 10), this test can be used to 
rule in favour of a particular seed-based model or, if N is found too large to 
rule out such models. For example, the number of textures predicted to have 
unwound in 10 x 10 degree MBR sky map between the time of recombination 
and today is less than 8. In fact, if reionisation is realized, as required by the 
texture model (Turok & Spergel 1990), this number will be much less than 
8. This implies that the proposed test may be efficiently used in this case 
since the predicted value of the kurtosis can be smaller by more than 10% 
compared to the Gaussian for 7 ~ I. 

Expressions similar to (35) may be easily obtained for higher moments 
of q. Using such expressions the proposed test can be applied even in cases 
where the signal to noise ratio 7 is not known, by using the measured con- 
straints on higher moments of q to eliminate 7. 
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It is straightforward to generalize the above analysis to higher dimensional 
cases. In fact we will show that the form of the generating function is the same 
in any number of dimensions. We will demonstrate the three dimensional 
case, applicable to large scale structure considerations. The two dimensional 
case corresponding to the MBR follows trivially from the three dimensional 
analysis. 

Consider a three dimensional rectangular area with coordinates x = 
(xi,X2,xs) such that —l<Xi<l{i = 1,2,3). In this case the wavenum- 

— * 

ber k becomes k = (hi, k 2 , ^3) and using the same analysis as in the one 
dimensional case it can be shown that 

P( Pl ,p 2 ) = R(pi,p 2 ) = ((^'/^^l^^l^^ae^^^^^^^y 

(36) 

Define now 



6 = 


7T 

jifaxi + k 2 x 2 + k 3 x 3 ) 


(37) 


6 = 




(38) 


& = 


jh%3 


(39) 



A change of variables from (xi,x 2 ,x 3 ) to (^1,^2,^3) leads to 

1 o 1 r+k 3 K r+k 2 Tr r+knT+^2+^3 ... f , . . f . 

P(pi,P2) = (— ) 3 TTT- / ^3/ d£ 2 ^ ie <tico S$1 +t 2 sm 6 ) 

111 k\k 2 k 3 J-k 3 n J-k 2 n J -kin+&+£3 

(40) 

which leads to a result identical to the one dimensional result (19) since 
ki,k 2 , k 3 are integers. It is trivial to see that the same is true for the two 
dimensional case. 

4 Discussion-The Power Spectrum 

So far we have investigated the statistical properties of the random func- 
tion Q(k) = qi(k) + iq 2 (k) which is only part of the Fourier modes of the 
perturbations. In fact we are interested in the full mode function 

g{k) = 9l {k)Q{k) (41) 
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Since the only random part of g(k) is Q(k), the statistical properties of g(k) 
are fully specified once we know such properties for Q(k). For example, the 
power spectrum of perturbations defined as 

P(k) =< \g(k)\ 2 > (42) 

is easily found using (29) to be 

P{k) = \g 1 {k)\ 2 <a\k)>=N\g 1 {k)\ 2 (43) 

for k ^ (obviously P(0) = N 2 |#i(£;)| 2 since Q(0) = N). It is possible to 
obtain the same result for the power spectrum by simply Fourier transforming 
the two point correlation function in coordinate space (Rice 1944; Scherrer 
& Bertschinger 1991). 

In conclusion, we have proposed a new statistical test for the identifica- 
tion of signatures of seed-based models in cosmological observations. The 
main advantages of investigating the statistics of perturbations in Fourier 
space rather than in coordinate space is that in Fourier space the analysis is 
valid for any geometry of superimposed seeds and can be directly applied to 
any particular experiment by simply selecting the Fourier modes for which 
the resolution and sky coverage of the experiment is most sensitive. No 
smoothing is needed as would be the case for coordinate space statistics. 

The statistical properties of the seed-like perturbations were shown to 
approach the Gaussian for a large number N of superimposed seeds. Thus, 
these properties can only distinguish efficiently models where the pertur- 
bations are produced by a small number of seeds per horizon scale. For 
an alternative statistic which can efficiently distiguish particular seed based 
models for larger N see Moessner, Perivolaropoulos & Brandenberger (1993). 

So far we have considered superposition of identical (but of any shape) 
seeds. However, our results can be easily generalized to variable seed magni- 
tude and extend in space. Such generalizations are shown in the Appendix. 
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Appendix A 



Generalizations 

An interesting generalization of our results can be provided by considering 
seeds of variable size. For example the size of perturbations induced by 
topological defects increases with cosmic time due to the growth of the hori- 
zon. Such effect can be taken into account by generalizing the Fourier space 
variable g 1 {k)Q{k) to 

M 

^ gi (Vk)Q(Vk) (44) 

j=0 

which corresponds to repeating the superposition of N seeds, M times while 
each time modifying the spatial scale of each seed by a factor of 2 in order 
to take into account the horizon growth (Vachaspati 1992; Perivolaropoulos 
1993b; Moessner, Perivolaropoulos & Brandenberger 1993). 

In this case the moment generating function P SU m{Pi,P2) for the sum 
of random variables is given (Feller 1971) by the product of the individual 
generating functions i.e. 

M 

P S um(pi,P2) = U(Jo(\gi(2 j k)\t)) N (45) 

3=0 

where the factor \g\{2^k)\ appears because we are now interested in the dis- 
tribution of the variable (44) as opposed to simply the variable Q(k). 

Finally, it is also straightforward to generalize our analysis to the case of 
seeds of variable magnitude. Such generalization would be needed in order 
to take into account the variable velocities of long cosmic strings. Consider 
for example the superposition of N seeds with Fourier transforms Xigi(k) 
(i = 1,...,N) where the coefficients Aj represent the relative magnitude of 
seed fluctuations. In this case the Fourier mode k becomes 

N 

g 1 (k)Y / ^e t ^=g 1 (k)Q x (k) (46) 

n=l 

and the generating function for the variable Q\{k) is 

Px(Pi,P2) = UM^) (47) 
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The above discussion is an attempt to show that our results are fairly gen- 
eral and can be easily adapted to the cases of particular seed-based models. 
Clearly, further work is needed to adapt the above analysis to any particu- 
lar model. Work in that direction for the cosmic string case is currently in 
progress. 



Figure Captions 

Figure 1: A comparison of the Gaussian (continous line) with P(q) for 
N = 10. 
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